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We investigate the chemical equilibration of the parton 
distributions in collisions of two heavy nuclei. We use initial 
conditions obtained from a self-screened parton cascade cal- 
culation and, for comparison, from the HIJING model. We 
consider a one-dimensional, as well as a three-dimensional ex- 
pansion of the parton plasma and find that the onset of the 
transverse expansion impedes the chemical equilibration. At 
energies of 100 GeV/nucleon, the results for one-dimensional 
and three-dimensional expansion are quite similar except at 
large values of the transverse radius. At energies of several 
TeV/nucleon, the plasma initially approaches chemical equi- 
librium, but then is driven away from it, when the transverse 
velocity gradients develop. We find that the total parton mul- 
tiplicity density remains essentially unaffected by the flow, 
but the individual concentrations of quarks, antiquarks, and 
gluons are sensitive to the transverse flow. The consequences 
of the flow are also discernible in the transverse momenta of 
the partons and in the lepton pair spectra, where the flow 
causes a violation of the so-called Mr scaling. 



I. INTRODUCTION 

Relativistic heavy ion collisions are being studied with 
the intention of investigating the properties of ultradense, 
strongly interacting matter. A central theme of this area 
of research is the search for a quark-gluon plasma (QGP) 
jjj . Extensive theoretical studies done over the last sev- 
eral years for nuclear reactions at collider energies have 
helped develop a picture of the colliding nuclei as two 
clouds of valence and sea partons which interact passing 
through each other ||. Partonic final-state interactions 
are then thought to produce a dense plasma of quarks 
and gluons. This plasma expands and becomes cooler 
and more dilute. If quantum chromodynamics admits a 
first-order deconfinement or chiral phase transition, this 
plasma will pass through a mixed phase of quarks, gluons, 
and hadrons, before the hadrons lose thermal contact and 
stream freely towards the detectors. 

What is the structure of the matter formed in these nu- 
clear collisions? Does the initial partonic system attain 
kinetic equilibrium? Probably yes, as the initial parton 
density is large, forcing the partons to suffer many col- 
lisions in a very short time ||. Does it attain chemical 
equilibration? This will depend on the time available 
to the partonic system before it converts into a 



mixed phase, or before perturbative QCD is no longer 
applicable. The time available for equilibration is per- 
haps too short (3-5 fm/c) at the energies (y/s < 100 
GeV/nucleon) accessible to the Relativistic Heavy Ion 
Collider (RHIC). At the energies (y/s < 3 TeV/nucleon) 
that will be achieved at the CERN Large Hadron Collider 
(LHC) this time could be large (more than 10 fm/c). If 
one considered only a longitudinal expansion of the sys- 
tem, the QGP formed at LHC energies could approach 
chemical equilibrium very closely, due to the higher ini- 
tial temperature predicted to be attained there. How- 
ever, then the life-time of the plasma would also be large 
enough to allow a rarefaction wave from the surface of the 
plasma to propagate to the center. As we shall discuss 
below, the large transverse velocity gradients developing 
in this process at LHC energies may drive the system 
away from the chemical equilibrium Q . 

It is evident that the evolution of the partonic system 
will depend on the initial conditions and the total time 
available to the system where the partonic picture still 
holds. In order to study this aspect, we use initial condi- 
tions from two models: the Self-Screened Parton Cascade 
(SSPC) model Jio) and the HIJING model jjjj. 

By now, a considerable progress has been achieved 
in our understanding of the parton cascades which de- 
velop in the wake of the collisions. Early calculations g] 
were done by assuming fixed px- and virtuality cut-offs 
for the partonic interactions to ensure the applicability 
of the perturbative expansion for QCD scattering pro- 
cesses. In the recently formulated self-screened parton 
cascade model early hard scatterings produce a medium 
which screens the longer ranged color fields associated 
with softer interactions. When two heavy nuclei col- 
lide at sufficiently high energy, the screening occurs on a 
length scale where perturbative QCD still applies. This 
approach yields predictions for the initial conditions of 
the forming QGP without the need for any ad-hoc mo- 
mentum and virtuality cut-off parameters. These calcu- 
lations also show that the QGP likely to be formed in 
such collisions could be very hot and initially far from 
chemical equilibrium. 

It is still useful to explore the dependence of the 
course of evolution of the plasma on the initial condi- 
tions. For this purpose we use the results of the HIJING 
Monte Carlo model which combines multiple hard or 
semihard parton scatterings with initial- and final-state 



1 



radiation with Lund string phenomenology for the ac- 
companying soft nonperturbative interactions to describe 
nuclear interactions. The uncertainties of the model are 
parametrized by allowing for variations of the initial fu- 
gacities. This approach has, for example, been used re- 
cently [l^ l in connection with the suppression of J/ip in 
an equilibrating partonic plasma. 

In the next section we briefly describe the hydrody- 
namic and chemical evolution of the plasma in a (1+1) 
dimensional longitudinal expansion and a (3+1) dimen- 
sional transverse expansion. The thermal photon and 
thermal lepton pair production is discussed in Section 
III, followed by a summary in Section IV. 



II. HYDRODYNAMIC EXPANSION AND 
CHEMICAL EQUILIBRATION 

A. Basic Equations 

We start with the assumption that the system achieves 
kinetic equilibrium by the time r^. We define the be- 
ginning of this epoch as the instant beyond which the 
momenta of the partons become locally isotropic. Local 
isotropy here is defined as the coincidence of the exponen- 
tial slopes of the longitudinal and transverse momentum 
distributions, when they are integrated over a comov- 
ing volume with linear dimensions of one partonic mean 
free path A/. This occurs after a proper time of about 
0.7 Xf ||, corresponding to t% ~ 0.2 — 0.3 fm/c under 
the conditions expected at collider energies Beyond 
this point, further expansion can be described by hydro- 
dynamic equations. The approach to the chemical equi- 
libration is then governed by a set of master equations 
which are driven by the two-body reactions (gg <-+ qq) 
and gluon multiplication and its inverse process, gluon 
fusion (gg <-» ggg). The hot matter continues to expand 
and cools due to expansion and chemical equilibration. 
Once the energy density reaches a critical value (here 
taken as ef = 1.45 GeV/fm 3 jlj|) we terminate the evo- 
lution before entering into the hadronization phase. 

Thus the expansion of the system is described by the 
equation for conservation of energy and momentum of an 
ideal fluid: 
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= A 
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T» v = (e + P)u»u v + Pg* v , (1) 



where e is the energy density and P is the pressure mea- 
sured in the frame comoving with the fluid. The four- 
velocity vector of the fluid satisfies the constraint 
u 2 = — 1. For a partially equilibrated plasma of mass- 
less particles, the equation of state can be written as Q 

e = 3P = [a 2 \ g + b 2 (X q + A 5 )] T 4 , (2) 

where a 2 = 8tt 2 /15, b 2 = 7n 2 N f /iO, N f ps 2.5 is the 
number of dynamical quark flavors, and A& is the fugac- 
ity for the parton species k. Here we have defined the 
fugacities through the relations, 



where is the equilibrium density for the parton species 
k: 
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C(3)T 3 = aiT 3 



n q = ^((3)N f T 3 = hT 3 . 



(4) 
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We further assume that X q = Xq. The equation of state 
(H) implies the speed of sound c s = 1/ 

We solve the hydrodynamic equations ([!]) with the as- 
sumption that the system undergoes a boost invariant 
longitudinal expansion along the z-axis and a cylindri- 
cally symmetric transverse expansion |l4|| . It is then suf- 
ficient to solve the problem for z — 0, because of the 
assumption of boost invariance. 

The master equations for the dominant chemical 
reactions gg <-> ggg and gg <-> qq are 

d^tigu^) = n g (R 2 ^3 - i? 3 ^2) - {rigRg^q - TlqRq^g) , 
d^n^) = d^nqu* 1 ) = n g R g ^ q - n q R q ^ g , (6) 

in an obvious notation. 

If we assume the system to undergo a purely longitu- 
dinal boost invariant expansion, (|l|) reduces to the well 
known relation Itt5fl 
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where r is the proper time. This equation implies 

e t 4/3 = const. (8) 
and the chemical master equations reduce to J3J 

= R 3 (1-X g )-2R 2 (l-^ 
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which are then solved numerically for the fugacities. The 
rate constants R 2 and -R3 are related to the rates appear- 
ing in (||) and are given by Q 

R 2 fa 0.24iV/c^A ff Tln(1.65/a s A ff ) ) 
R^l.2alT{2Xg-X 2 g ) 1 /\ (10) 

where the color Debye screening and the Landau- 
Pomeranchuk-Migdal effect suppressing the induced 
gluon radiation have been taken into account, explicitly. 

In case of transverse expansion, the master equations 
can be shown p| to lead to partial differential equations: 
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where v r is the transverse velocity and 7 = 1/y 1 — tc- 
It is easy to verify that the equations ( Jll| ) reduce to (g) 
in the absence of transverse expansion. 

The hydrodynamic equations (Q) are solved numeri- 
cally to get e(r, t) and i> r (r, t), which serve as input into 
the equations (^Tj) for the fugacities. In all our solutions 
we have assumed that the initial transverse velocities are 
zero. We have verified that our results near r = closely 
follow the results for a purely longitudinal expansion, till 
the time when the fluid is disturbed by the rarefaction 
wave traveling from the surface to the center. 



B. Results for RHIC Energies 

In the following, we shall discuss the results obtained 
for the initial conditions listed in Table I, which are pre- 
dicted for a central collision of gold nuclei. We shall 
discuss the results for the initial conditions from SSPC 
calculations in some detail, and then present the final 
results for the initial conditions from the HIJING model. 

Let us first look at conditions likely to prevail at RHIC 
energies. For the case of a boost-invariant and a purely 
longitudinal expansion, the results are discussed in the 
literature Jig ]. Here we concentrate on the question, how 
the transverse expansion of the plasma affects the chemi- 
cal equilibration, and how the consequences of the trans- 
verse flow can best be identified. 

Recall that for a boost-invariant longitudinal expan- 
sion, Eq.(|l]) provides that er 4 / 3 is a constant (||). Let us 
look at the constant energy density contours, 



e(r,t) = e 4 /iV 4/3 



(12) 



for N — 1,2,3, ...,16 (Figure la). The choice of hy- 
persurfaces of equal comoving energy density is natural, 
because it interpolates smoothly between the initial con- 
dition and the final hadronization hypersurface, where 
the quark-gluon plasma freezes out into hadrons. If there 
were no transverse expansion, all these contours would be 
parallel to the line for N = 1 and given by 



extending up to r = Rt- 

For the initial energy density assumed at RHIC, 
N = 16 corresponds to an energy density of about 1.5 
GeV/fm 3 , which coincides with our choice of the density 
tf where the plasma finally hadronizes. We have also 
shown a line r = Rt — c s t which indicates the radial size 
of the region not yet affected by the rarefaction wave at 
any given time t. We see that, as expected, the fluid 
beyond r = 4 fm is likely to be affected by the flow, by 
the time the system has cooled to the edge of the QGP 
phase. 

In Fig. lb-d we have given the gluon fugacity, the 
quark fugacity, and the transverse velocity along these 
contours. We see that deviations from the expectations 
for a longitudinal expansion [ fl6| , which we closely repro- 
duce for small r, have a clear origin in the growth of the 
transverse velocity. 
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FIG. 1. (a) Constant energy density contours for a trans- 
versely expanding quark-gluon plasma likely to be created in 
Au+Au collisions at RHIC. For a boost invariant longitudinal 
expansion these contours would be equidistant lines parallel 
to N = 1. The dashed line gives the distance from the axis 
where the rarefaction wave has arrived at time t, and beyond 
which the fluid is strongly affected by the flow, (b) Gluon fu- 
gacity along the constant energy density contours, (c) Quark 
fugacities along the constant energy density contours, (d) 
Transverse velocities along the constant energy density con- 
tours. 

For small values of N, we see that the plasma at the 
center has more time at its disposal to develop while 
at large values of r the time available to the plasma is 
shorter, leading to a decrease in the fugacities (at con- 
stant e) as r increases. However, this is not the most 
interesting observation. We see a unique feature, which 
is quite evident at large r but is beginning to emerge 
at smaller r as well, especially for the gluon fugacity: 
The fugacities initially increase with TV, but then start 
decreasing. Thus the region where transverse velocities 
are large, the plasma ends its journey further away from 
chemical equilibrium, as compared to the region where 
the transverse velocity is still small. We shall see later 
that this has its origin in the velocity gradients which 
develop in the fluid. 

This has an interesting consequence, namely, the 



plasma nearer to the surface will be hotter than in the 
interior, for a given energy density. If we make the rea- 
sonable assumption that the transition point between the 
quark-gluon plasma and the hadron phase is determined 
by a certain critical value of the energy density, then 
our result implies that the outer regions of the cylindri- 
cal reaction zone hadronize at a higher temperature than 
the inner region. In other words, it is quite likely that 
the hadronic matter may be formed at different temper- 
atures in different regions. This enhances the complexity 
of the final state produced in the hydrodynamic model 
compared to what was expected earlier. 

How does the flow affect the momentum distribution 
of the partons? We can evaluate the evolution of the 
parton distributions with time for the longitudinal ex- 
pansion, and compare it with the corresponding results 
for the transverse expansion, using the Cooper-Frye for- 
mula 



dN 



fJ 



dN 



d 3 p/E d 2 p T dy (2tt) 



f{x,p)p^da^ 



(14) 



where f(x,p) is the phase-space distribution of the par- 
ton, g is the color, spin, and flavor degeneracy, and a is 
the surface described by the contours given by ( [l^ ) or 
(|l3|). For massless particles and in the absence of trans- 
verse flow simplifications arise, which we shall not discuss 
here Q. 

The resulting transverse momentum spectra for glu- 
ons and quarks is shown in Figs. 2a~b. We see that the 
parton distributions at large px are affected strongly as 
the flow develops. It is known that if the plasma were to 
undergo a first order phase transition, then the nonvan- 
ishing transverse velocity at the beginning of the mixed 
phase will reduce its life-time. We shall see, however, that 
due to several other competing factors, like lower tem- 
perature at later times, the flow only marginally affects 
the production of thermal photons and lepton pairs from 
the quark-gluon phase at RHIC energies. The presence 
of significant collective transverse flow is clearly visible 
in the positive curvature of the p-r-spectra at late times 
(N = 12,16). The slope of these spectra at large pr, 
resulting from a superposition of thermal and collective 
motion, agrees remarkably well with the homogeneous 
slope of the initial spectrum (N = 1) and thus serves 
as an indicator of the very high apparent temperature 
achieved at the moment of thcrmalization. 
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FIG. 2. Evolution of the momentum distributions of (a) 
gluons and (b) quarks, and (c) of the multiplicity of partons 
with (solid curves) and without (dashed curves) transverse 
flow at RHIC. 



quarks (and antiquarks) whereas the number of gluons 
decreases. 

However, the detailed composition of the partonic mat- 
ter at the end of the QGP phase is different for the 
case with flow. We have approximately 12% less quarks 
and antiquarks but about 10% more gluons compared 
to our result without flow. It is difficult to predict how 
these rather subtle differences in the composition of the 
plasma will affect the final particle yields. We may spec- 
ulate that an increased gluon abundance could result in 
the enhanced production of hadrons with a significant 
"valence-glue" component, such as r\- and 7y'-mesons or 
even glueballs. 



C. Results for LHC Energies 

Let us next turn to the results for Au+Au collisions 
at LHC energies, some of which have been reported ear- 
lier j^]. In Figs. 3a-d, we present the constant energy 
density contours, the fugacities for gluons and quarks, 
and the transverse velocities. The representation is the 
same as the one for the RHIC-related figures, except that 
we have plotted only a reduced number of equal-energy 
contours. Due to the very large number of contours be- 
tween the initial energy density and the hadronization 
density e/ = 1.45 GeV/fm 3 we have plotted only the 
contours TV = 1, 6, 11, . . . , 71 in the energy density plot 
and TV = 1, 11, 21, . . . , 71 in the plots for the other three 
quantities. 
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We can also evaluate the multiplicity density of the 
partons as a function of N to see the increase in entropy, 
as the chemical equilibration proceeds, and investigate 
its dependence on the flow. We see from Fig. 2c that the 
total entropy at the end of the QGP phase is essentially 
identical for the scenarios with and without transverse 
flow. We also see that as the chemical equilibration is 
accompanied by a steady increase in the net number of 
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FIG. 3. (a)-(d): Same as Fig. 1, at LHC. 

The same effects pertaining to the collective transverse 
flow are observed as at the lower RHIC energy, but they 
are more pronounced. Due to the longer lifetime of the 
plasma phase at LHC energies, the effects of the trans- 
verse expansion are generally more dramatic. The rar- 
efaction wave reaches the center after about 12 fm/c, and 
all of the matter participates in the transverse flow at the 
moment of hadronization. The more rapid cooling of the 
system is also indicated by the decreasing separation of 
the contours at later times. This affects the chemical 
equilibration process in several competing ways. 



We note that, as before, the quark fugacities lag be- 
hind the gluon fugacities at all times and all radial dis- 
tances. This is not surprising, considering the lower start- 
ing value for X q . We also see again that the fugacities in 
the interior, at first, increase with passage of time. How- 
ever, beyond N = 31, corresponding to about 8 fm/c near 
the center, when the transverse flow has been substan- 
tially developed, the fugacities begin decreasing again, 
except at r = 0, where they continue to increase, as all 
the radial derivatives in (^Tj) are small there. 

Let us pause to understand this unexpected feature. 
According to (^), the evolution of the parton densities is 
governed by d^{nu^) which consists of two terms: u^d^n 
gives the rate of change of the density in the comoving 
frame, and nd^u^ describes the rate of change of the 
density due to the expansion of the fluid clement in the 
comoving frame JlSf] . Once the transverse expansion of 
the fluid starts developing, the second term and also the 
radial derivatives in ( |TT| ) grow very rapidly and drive the 
system away from chemical equilibrium. 

We thus see that two mechanisms potentially con- 
tribute to the depletion of the parton densities: expan- 
sion and flow. Let us concentrate on the consequences 
of the latter, as the consequences of the expansion alone 
are easy to visualize. Figure 4 depicts a fluid element 
at a distance r from the axis at some instant of time, in 
a central collision for the conditions of: (a) no flow, (b) 
uniform transverse flow, and (c) a differential transverse 
flow with the velocity field v r (r, t) as obtained in the 
present work. When there is no flow, the thermal motion 
of the partons brings the same number of partons into the 
volume element as are leaving it on all sides. Thus any 
production of new partons due to the chemical reactions 
in the volume element will increase the fugacities, which 
will hence rise monotonically since the right hand sides 
of the master equations (0,fLl]) remain positive [EUlG] . 
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FIG. 4. Depletion of partons for, (a) no flow, (b) uniform 
flow, and (c) flow with velocity gradient in a fluid element at 
r. 

Next consider the condition of a uniform transverse 
flow, which is sometimes used as a simplifying assump- 
tion |2(J. This scenario still does not produce a net out- 
flow of partons from the volume element, because the loss 
of partons on one side is exactly balanced by the gain of 
partons flowing in from the opposite side. Thus, if there 
is a creation of new partons due to chemical reactions 
then, as before, the fugacities will increase with time. 

However, when the transverse velocity is driven by the 
pressure gradient as in our calculation, v r (r+dr) > v r (r). 
Then the velocity gradient will ensure that more partons 
are leaving the fluid element at r + dr than are entering 
at r. This leads to an additional depletion of the par- 
ton density, beyond the dilution caused by the overall 
expansion. If the velocity gradient is sufficiently large, 
the production of new partons by the chemical reactions 
cannot cope with the depletion and the fugacities will 
start to fall. The richness of this scenario is completely 
absent when the simplifying assumptions of no-flow or of 
a uniform transverse flow are made. 

Let us return to the discussion of the consequences of 
the transverse flow on the momentum distribution of the 
partons at LHC energies. In Figs. 5a-b, we show the evo- 
lution of the gluon and quark momentum distributions 
over the course of the expansion. We see that the distri- 
bution at later times is strongly affected by the flow. It 
is again apparent that there is an approximate balance 
between the decreasing temperature and the increasing 
transverse velocity which keeps the parton distributions 
at high pt largely independent of TV. This is reminis- 
cent of a similar behavior observed in pion spectra when 
transverse flow is introduced [|TJ . 
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FIG. 5. (a)-(c) Same as Fig. 2, for LHC. 

Figure 5c shows our results for the evolution of the 
parton multiplicity. We see a significant production of 
quarks and antiquarks, in fact an increase by a factor 
of about 2.5, and a decrease in the number of gluons by 
about 25%, as a result of the chemical reactions in the 
plasma. This results in a net increase of the parton mul- 
tiplicity by approximately 40%. As before at RHIC en- 
ergies, we find that the transverse flow does not strongly 
affect the net production of partons, although it influ- 
ences the individual abundances of gluons and quarks. 

It is probably difficult to probe the off-chemical equilib- 
rium nature of the light quark and gluon fractions in the 
plasma. However, we can think of a practical way of test- 
ing the chemical equilibrium in the cases of heavier quark 
flavors, viz., strange, charm, and even bottom quarks, 
by inferring this information from the ratios of the (va- 
lence) quark content of the emitted hadrons (2^j23|] . We 
can analogously define the ratio of heavier quarks to light 
quarks as 
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where Q stands for one of the heavy quarks. It is not 
unlikely that this ratio can be largely determined from 
the final production of strange, charm, or bottom mesons 
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and pions 0. 

Simulations of the pure SU(3) lattice gauge theory have 
shown that the speed of sound drops significantly as the 
temperature approaches the phase transition tempera- 
ture from above p4| |. If this observation carries over 
to full QCD, our results could be modified. A lower 
speed of sound generally leads to a reduced production 
of transverse flow and, hence, to an increased life-time of 
the QGP. This effect should be especially pronounced at 
RHIC energies pq] . 

Can we devise means for identifying the large trans- 
verse flow velocities predicted here? As mentioned ear- 
lier, the life-time of the mixed phase, which follows the 
QGP phase, gets shorter if the transverse velocity in- 
creases. Because the transverse flow will remain frozen 
in the mixed phase, the final hadronic flow pattern should 
be a rather reliable measure of the transverse velocity es- 
tablished during the QGP phase. In a forthcoming publi- 
cation, we shall report on the effect of the transverse flow 
on the pT-spectra of charmed hadrons which are consid- 
ered to be good probes of flow due to their large mass. In 
the next section we look at the spectra of single photons 
and lepton pairs and see how far could they be affected 
by these consideration. 



D. Other Initial Conditions 

Let us now briefly examine the consequences of us- 
ing different initial thermalization conditions. We 
only present the results for the hadronization contour, 
e(r, t) = 1.45 GeV/fm 3 , and for the gluon and quark fu- 
gacities attained along this contour for RHIC (Fig. 6a-b) 
and LHC (Fig. 7a-b) energies. We select several initial 
conditions obtained by modifying predictions of the HI- 
JING model, as listed in Table I. The original prediction 
of the HIJING model jll] is denoted as I. The second 
set (II) of initial conditions is obtained by increasing the 
initial fugacities, and hence the energy density, by a fac- 
tor four. The third set (III) is obtained by retaining the 
fugacities of the second set, but decreasing the temper- 
ature to account for the production of soft partons from 
the color field, by some arbitrary amount. It is of inter- 
est to recall that the Parton Cascade Model of Geiger || 
predicts an energy density of about 120 GeV/fm 3 , and 
a temperature of about 590 MeV, at t = 0.25 fm/c in 
central collisions of two gold nuclei at RHIC energies. 
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FIG. 6. (a) Constant energy density contours for e = 1.45 
GeV/fm 3 for the initial conditions obtained from HIJING, 
given in Table I. (b) Gluon and quark fugacities along these 
contours. 
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FIG. 7. Same as Fig. 6, for LHC. 

If we use the original predictions (I) of the HIJING 
model, then the life-time of the QGP phase is very small 
(less than 2 fm/c), and the matter at the moment of 
hadronization is extremely dilute compared with chemi- 
cally equilibrated matter. Changing the initial fugacities 
by a factor of 4 increases the life-time of the QGP phase 
considerably and also brings the system close to chemi- 
cal equilibrium. Note that the initial fugacities differ by 
one order of magnitude between the conditions I and III 
at RHIC energies. Note also that the final fugacities at 
large r are quite different, even though the energy den- 
sities are identical, signaling very different temperatures 
at the end of the QGP phase. 

Results for LHC energies indicate the effect of trans- 
verse expansion for r > 4 fm, for the original predictions 
of HIJING, marked as I, and the entire fluid is affected 
by the flow for the other initial conditions, which assume 
much larger initial energy densities and fugacities. The 
final fugacities for the conditions II and III are essen- 
tially identical and show the effect of evolution away from 
chemical equilibrium, as the flow velocity gets large. The 
transverse dimension of the system grows by almost 1 fm 
during the QGP phase itself for the initial conditions II 
and III. Since the fugacities shown in these two figures 
are for the same energy density e / , we note that the final 
temperature for the standard case (I) will be much larger 
than in the scenarios II and III. 



III. THERMAL PHOTONS AND LEPTON PAIRS 

A. Photon Spectra 

Thermal photons and lepton pairs are primary probes 
of all the stages of the nuclear reaction. They are also 
expected to carry valuable imprints of the transverse ex- 
pansion of the system. Thermal photons from the quark- 
gluon plasma have their origin in the Compton (qg — ► (77) 
and annihilation processes (qq — > gj). It is useful to re- 
call the expressions for the rates for these contributions 



for a non-equilibrated quark-gluon plasma pq|. For the 
approximations of the fugacities that we have used, the 
Compton process contributes with the rate 
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and the rate of the radiative annihilation process is, 
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Here C = 0.577721. . . , e q is the electric charge of the 
quark and the parameter k c is related to the thermal 
mass of the quarks in the medium. 

Recalling that our plasma is gluon-rich and quark- 
poor, we anticipate that the yield of the thermal photons 
would be dominated by the Compton contribution. In 
Figs. 8 and 9, we give our results for RHIC and LHC 
energies, respectively, for the initial conditions obtained 
from the SSPC fllC] . We show the Compton and annihi- 
lation contributions separately for the scenario involving 
transverse expansion. The results for longitudinal expan- 
sion are given for comparison. 
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FIG. 8. Distribution of thermal photons from the QGP 
phase at RHIC. The Compton and annihilation yields, and 
their sum are shown for the case with transverse expansion. 
Results are also given for a purely longitudinal flow. Prompt 
photons, whose production is governed by structure functions, 
are seen to dominate the yield for pr > 3-4 GeV. 
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FIG. 9. Same as Fig. 8 at LHC. The prompt photons from 
fragmentation of quark jets are also shown. 

We see that the transverse expansion of the fluid leaves 
the thermal photon production from the QGP at RHIC 
energies essentially unchanged for pr > 1 GeV. The re- 
duction of the thermal photon production at lower px will 
be difficult to identify, as it will be overshadowed by pho- 
tons from hadronic reactions and hadronic decays. We 
also see that prompt photons [^7j will dominate the yield 
beyond pr > 4 GeV. At first this result looks surprising, 
in view of the large modifications of the final parton spec- 
tra seen earlier (cf. Figs. 2 and 5L However, note that 
the rates for photon emission ( |l6| , [l7| ) carry an additional 
weight factor of T 2 , which suppresses the contributions 
from later times. Note also that the time available for 
radiation at a given temperature is significantly reduced 
by the transverse expansion. Overall, we conclude that 
the thermal photon production from the plasma phase 
is remarkably insensitive to the transverse expansion at 
RHIC energies. 

The production of thermal photons with large px re- 
mains essentially unchanged at LHC, as well, as they 
have origin in the early hot stages, when the flow effects 
are still small. On the other hand, the production of low 
Pt photons decreases considerably, due to reduction in 
the space-time volume occupied by colder matter in the 
presence of transverse expansion. While the prompt pho- 
ton production p7[ may remain lower than the thermal 
photon yield up to px — 5 GeV, the background contri- 
bution of photons fragmented off high-p<r quark jets is 
large. 

B. Lepton Pair Spectra 

The mass distribution of lepton pairs obtained from 
quark annihilation are given in Figs. 10 and 11, for RHIC 
and LHC energies, respectively, without and with trans- 
verse expansion. The Drell-Yan contribution |^8| is also 
shown. We again see a negligible effect of the flow on the 
mass distribution of lepton pairs at RHIC energies, but a 



factor of 2-3 reduction is seen in the yield at LHC ener- 
gies at invariant masses of 1-2 GeV, due to the reduced 
life-time of the plasma. We also note that the Drell-Yan 
contribution dominates the yield beyond M > 4 GeV 
at LHC, whereas this zone shifts to much smaller M at 
RHIC. This is primarily due to very low quark densities 
of the plasma. 
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FIG. 10. Mass distribution of thermal lepton pairs from 
the QGP phase at RHIC with (solid curves) and without 
(dashed curves) transverse flow. The Drell- Yan contribution 
is also given. 
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FIG. 11. Same as Fig. 10 for LHC. 

These small deviations in the production of thermal 
photon and lepton pairs due to the transverse expansion 
are rather astounding. We have already noted that this 
has its origin in the strong temperature dependence of the 
production rates, combined with the accelerated cooling 
of a transversely expanding plasma. This raises the ques- 
tion whether it is at all possible to identify the presence 
of transverse flow from the lepton pair spectra. Recall 
that the transverse mass distribution of the lepton pairs 
from the QGP in presence of a flow will be given by, 
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In the absence of a flow, this reduces to 
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which scales with Mt in a characteristic manner which 
is known as Mt- scaling ]2!||3(|. This scaling will be 
violated in the presence of a transverse flow ||30| . 

In Figs. 12 and 13 we present our predictions for the 
ratios of the transverse mass distributions for M = 1 
GeV to that for M = 2, 3, and 4 GeV at RHIC and 
LHC energies, respectively. As indicated earlier these 
ratios are identically one for a purely longitudinal expan- 
sion. The deviations predicted here are not large, but we 
hope that they will be observable in the high statistics 
data that will be forthcoming from RHIC and LHC. Ad- 
ditional information about the parton fugacities is con- 
tained in the absolute yields of the photon and lepton pair 
spectra, because the thermal photon yield is dominated 
by the Compton processes (oc X g X q ), whereas the lepton 
pair yield is mainly due to quark-antiquark annihilation 
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FIG. 12. Ratio of transverse mass distribution for thermal 
lepton pairs for M = 1 GeV to that for a given M. In absence 
of flow, it is identically unity. 
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IV. SUMMARY 

We have studied the evolution of an equilibrating and 
transversely expanding quark-gluon plasma, with ini- 
tial conditions that have been obtained either from the 
self-screened parton cascade model or from the HIJING 
model. We have compared the results from a boost in- 
variant longitudinal, as well as from a transverse expan- 
sion of the plasma. The velocity gradients generated due 
to flow were found to affect the individual concentrations 
of the quarks, antiquarks, and gluons, without affecting 
the total final entropy. 

Large transverse velocities are developing by the end 
of the QGP phase, especially at LHC energies. The flow 
is also seen to drive the system away from chemical equi- 
librium. As the process of chemical equilibration accel- 
erates cooling, and as the rates of emission of photons 
and lepton pairs depend strongly on the temperature, 
the final spectra are only slightly affected by them. This 
has its origin in the fact that photons and lepton pairs 
originate from the early hot stage when the flow is still 
small. However, the transverse flow is seen to introduce 
a slight violation of the so-called Mt scaling of the lep- 
ton pair spectra. Likely consequences on the post-QGP 
epoch could be a copious production of low pt pions from 
the fragmentation of residual gluons. We also note that 
the temperature of the matter at the end of the so-called 
QGP phase, may be higher at the surface than in the 
center of the reaction volume. 
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